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Abstract 

We present a calculation of corrections to the lifetime differences of Bg mesons AF^^ in the 

heavy-quark expansion. We find that they are small to significantly affect AF^^ and present the result 
for lifetime difference including non-perturbative l/m^, and l/m^ corrections. We also analyze the generic 
AB = 1 New Physics contributions to AF^^ and provide several examples. 
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I. INTRODUCTION 



Meson- antimeson mixing serves as an indispensable way of placing constraints on various models 
of New Physics (NP). This is usually ascribed to the fact that this process only occurs at the 
one-loop level in the Standard Model (SM) of electroweak interactions. This makes it sensitive 
to the effects of possible NP particles in the loops or even to new tree-level interactions that 
can possibly contribute to the flavor-changing AQ = 2 interactions. These interactions induce 
non-diagonal terms in the meson- antimeson mass matrix that describes the dynamics of those 
states. Diagonalizing this mass matrix gives two mass eigenstates that are superpositions of flavor 
eigenstates. In the Bg system mass eigenstates, denoted as "heavy" \Bh) and "hght" \Bl), 

\BH)^p\B,) + q\Bs), 

\BL)=p\Bs)-q\B,), (1) 

were predicted to have a rather significant mass and width differences, 

AMb, ^Mh- Ml, AFb, = - F^, (2) 

where M^^, and Th^l denote mass and lifetime differences of mass eigenstates. Since in the 
Standard Model the mass difference is dominated by the top quark contributions, it is computable 
with great accuracy. Thus one might expect that possible NP contributions can be easily isolated. 
Unfortunately, a recent observation of mass difference of mass eigenstates in Bg mixing by CDF [1] 
and DO [2], 

AMs^ = 17.77 ±0.10 ±0.07 ps"^ (CDF), 

17 ps-^ < AMb, < 21 ps-^ (DO), (3) 

put the hopes of spectacular NP effects in Bg system rest. In fact, analyses of mixing in the strange, 
charm and beauty quark systems all yielded positive signals, yet all of those signals seem to be 
explained quite well by the SM interactions. Yet, some contribution from New Physics particles 
is still possible, so even the energy scales above those directly accessible at the Tevatron or LHC 
can be probed with Bg mixing, provided that QCD sum rule [3] or lattice QCD [4-7] calculations 
supply the relevant hadronic parameters with sufficient accuracy. 
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In addition to the mass difference AM^^ , a number of experimental collaborations reported the 
observation of a hfetime difference AF^^ in the Bg system. Combining recent result from DO [8] 
with earher measurements from CDF [9] and ALEPH [10], Particle Data Group (PDG) quotes [11] 

AF^, = 0.161°:}° ps-\ ^ = 0.12ll°:°«3, (4) 

while Heavy Flavor Averaging Group (HFAG) [12] gives 

AF^, = O.OTllHfr PS-\ ^ = O.lOit',^. (5) 

Differently from the mass difference AMb^, the lifetime difference AF^^ is definitely dominated 
by the SM contributions, as it is generated by the on-shell intermediate states [13-16]. While 
this might appear to make it less exciting for indirect searches for New Physics, besides "merely" 
providing yet another test for heavy quark expansion, it is nonetheless a useful quantity for a 
combined analysis of possible NP contributions to Bg^ — -65° mixing [17-20]. 

It has been argued [19] that CP-violating NP contributions to AB = 2 amphtudes can only 
reduce the experimentally-observed lifetime difference compared to its SM value, therefore it is 
important to have an accurate theoretical evaluation of AF^^ in the SM. It is also important to 
note that AB — 1 NP contributions can affect AF^^, but do not have to follow the same pattern. 
Indeed, the level at which AB — 1 NP can affect AF^^ depends both on the particular extension of 
the SM, as well as on the projected accuracy of lattice calculations of hadronic parameters which 
drives the uncertainties on the theoretical prediction of AF^^. So it is advantageous to evaluate 
the effect of NP contributions. 

This paper is organized as follows. We set up the relevant formalism and argue for the need 
to compute 1/ml corrections to leading and next-to-leading effects in Sect. II. In Sect. Ill we 
discuss the impact of 1/ml corrections to the hfetime differences of Bg mesons and assess the 
convergence of the l/m^ expansion. We also present the complete SM results for AF^^ including 
1/ml corrections. We then discuss the possible effects from AB — 1 New Physics contributions in 
Sect. IV. Finally, we present our conclusions in Sect. V. 

II. FORMALISM 

In the limit of exact CP conservation the mass eigenstates of the B^-B^^ system are \Bh/l) — 
{\Bs) ± \Bs))/V^, with the convention CP\Bs) — — The width difference between mass 
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eigenstates is then given by [13] 

afb, = - r^^ = -2 = -2 , (6) 

where Tij are the elements of the decay-width matrix, i,j = 1,2 (|1) = |-Bs), |2) = \Bs)). 

We use the optical theorem to relate the off-diagonal elements of the decay-width matrix F 
entering the neutral S-meson oscillations to the imaginary part of the forward matrix element of 
the transition operator T: 

r2i(S,) = ^^{B,\r\B,), T^lmiJ d''xT{H,s{x)H,s{0)} . (7) 

Here Heff is the low energy effective weak Hamiltonian mediating bottom-quark decays. The 
component that is relevant for reads explicitly 



(8) 



defining the operators 

Ql = {biCj)v-A{CjSi)v-A, Q2 = {kCi)v-A{CjSj)v-A, (9) 

Q3 = {biSi)v-A{qjqj)v-A, Qa = (biSj)v-A{qjqi)v-A, (10) 

Qb = {biSi)v-A{qjqj)v+A, Qe = {biSj)v-A{qjqi)v+A, (n) 

Q8 = ^m,6,a^-(l-75)i;>,G^.. (12) 

Here i,j are color indices and a summation over q — u, d, s, c, b is implied. V ± A refers to 
7'^(1 ± 75) and S — P (which we need below) to (1 — 75). Ci, . . . , Ce are the corresponding Wilson 
coefficient functions at the renormalization scale which are known at next-to-leading order. We 
have also included the chromomagnetic operator Qs, contributing to T at 0(tts). Note that for a 
negative Cg, as conventionally used in the literature, the Feynman rule for the quark-gluon vertex 
is —ig^fj,T^. A detailed review and explicit expressions may be found in [21]. Cabibbo-suppressed 
channels have been neglected in Eq. (8). 

In the heavy-quark limit, the energy release supplied by the b-quark is large, so the correlator 
in Eq. (7) is dominated by short-distance physics [22]. An Operator Product Expansion (OPE) 
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can be constructed for Eq. (7), which results in its expansion as a series of matrix elements of local 
operators of increasing dimension suppressed by powers of l/m^: 

r2i(5.) = Y.{Bs\%\B,) = E ^{B,\Ot''=\^,)\B:). (13) 

In other words, the calculation of r2i(-Bs) is equivalent to computing the matching coefficients of 
the effective AB — 2 Lagrangian with subsequent computation of its matrix elements. Eventually 
the scale dependence of the Wilson coefficients in Eq. (13) is bound to match the scale dependence 
of the computed matrix elements. 

Expanding the operator product (7) for small x ~ i/rrib, the transition operator T can be 
written to leading order in the l/mj expansion as [13, 14] 

r = -^(K;KJ' [F{z)QM + Fs{z)Qs{f^,)] , (14) 

which results in [15] 

X { [(1 - z) (2 C1C2 + N,Cl) + (1 - ^z)Cf/2] (Q) + (1 + 2z) (2 C1C2 + N,Cl C^) {Qs)}(,l5) 

with z — ml /ml and the basis of AB — 2 operators-*^ 

Q = {biSi)v-A{bjSj)v-A, Qs = {biSi)s-p{bjSj)s-p ■ (16) 

In writing Eq. (14) we have used the Fierz identities and the equations of motion to eliminate the 
color re-arranged operators 

Q = {biSj)v-A{bjSi)v-A: Qs = {biSj)s-p{bjSi)s-p , (17) 

always working to leading order in l/m^. Note that (...) denote matrix elements of the above 
operators taken between Bg and Bs states. The Wilson coefficients F and Fs can be extracted by 
computing the matrix elements between quark states of T in Eq. (7). 

The coefficients in the transition operator (14) at next-to-leading order, still neglecting the 
penguin sector, can be written as [14]: 

F{z) = Fn{z)ClM + F,2{z)CMCM + F^2{z)Cl{^x,) , (18) 



^ It was recently argued that better-converging results can be obtained in a modified basis [16]. 
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(20) 
(21) 



F.,{z) = F^{z) + ^Fi^\z), (19) 
and similarly for Fs{z). The leading order functions F^'^\ F^J^ read explicitly 

F[f{z) = 3y/T^ (l-z), F^%{z) = 3^1^ (1 + ^z), 

f[^{z) = 2x/r^ {l-z), F^%{z) = 2VT^ (1 + 2z), 
Fi?i^) = ^(1 - ^zf', Ffl,{z) = -Vl^z (1 + 2z). (22) 

The next-to-leading order (NLO) expressions Flj \ Fg^Jj are given in Ref. [14]. 
The penguin correction to Eq. (14), 

% = -^iv:,Ksr [Piz)Q + Ps{z)Qs] , (23) 

is also shown in Ref. [14]. 



III. SUBLEADING 1/m^ CORRECTIONS 

Here we present the higher order corrections to r2i{Bs) in Eq. (15) in the heavy-quark expansion, 
denoted below as Si/m and 61/^2: 

(24) 

The matrix elements for Q and Qs are known to be [13-15] 

{Q)=(B,\Q\B,) = flMl2 (1 + i-) B, 



iQsMB.lQsm = -flM,^^^^ (2 - -) B,, (25) 

Si/m^{Bs\%./m\Bs), and (5i/^2 = (Ss|Tl/^2 

where Mb^ and /b^ are the mass and decay constant of the Bs meson and Nc is the number of colors. 
The parameters B and Bs are defined such that B — Bg — 1 corresponds to the factorization (or 
'vacuum insertion') approach, which can provide a first estimate. 
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FIG. 1: Calculation of kinetic l/m\, and corrections. The operators of Eqs. (26) 
and (30) are obtained by expanding the diagrams in powers of light quark momen- 
tum. 



A. l/rrib corrections 

The 1/mb corrections are computed, as in Ref. [13, 15, 22, 23], by expanding the forward 
scattering amplitude of Eq. (7) in the light-quark momentum and matching the result onto the 
operators containing derivative insertions (see Fig. 1). The contributions can be written in 
the following form: 

T^/m^sfl^z I (1 + 2z) [Cl (i?2 + 2i?4) - 2 (2 C1C2 + N^Cl) (i?i + R2) 



1-Az 

where the operators Ri are defined as 



(2 C1C2 + N,C^) (i?2 + 2R3) + 2 i?3 



(26) 



mt 





1 T 




—b. 




mi 






R4 — 


l-i 




nib 



Their matrix elements read [13, 15]: 

1 \ m 



{Bs\Ri\Bs)^{2 + 



{B,\R,\B,)^l + ^) /|M|J-^-l) S|, (S,|it:4|S,) = -/|M|J-f -1) 5|.(28) 



2Nr 



mt 



Some of those parameters have been computed in lattice QCD [4-7].' In this paper we use the 
results of Ref. [4]. 



^ For estimates of these matrix elements based on QCD sum rules, see Ref. [3] . 
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FIG. 2: 1/m^-corrections from gluonic operators. 



The color-rearranged operators Ri that follow from the expressions for Ri by interchanging color 
indexes of bi and Sj Dirac spinors have been eliminated using Fierz identities and the equations 
of motion as in Eq. (16). Note that the above result contains full QCD 6-fields, thus there is no 
immediate power counting available for these operators. The power counting becomes manifest at 
the level of the matrix elements. 

B. i/Tri^ corrections 

It was shown in Refs. [13, 15] that l/m(,-corrections are quite large, so it is important to assess 
the convergence of 1/mfe-expansion in the calculation of the Bg lifetime difference. In order to do so, 

we compute a set of 5i/m2 corrections to leading order. As expected, at this order more operators 
will contribute. We will parametrize the 1/ml corrections similarly to our parametrization of 
l/mb effects above and use the factorization approximation to assess their contributions to the Bs 
lifetime difference. 

Two classes of corrections arise at this order. One class involves kinetic corrections which can be 
computed in a way analogous to the previous case by expanding the forward scattering amplitudes 
in the powers of the light-quark momentum. A second class involves corrections arising from the 
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interaction with background gluon fields. The complete set of corrections is the sum of those, 



rr rrkin , 'TG 

Let us consider those classes of corrections in turn. The kinetic corrections can be written as 
T^t,^Vl^z\- —^(3 - lOz) \CIW^ + (2 C1C2 + N,Cl){W^ + W^/^) 



(29) 



12^2 rn^ 



1 — izml 



- (2 C1C2 + N,C',){Qs + g/2)] 



24^2 



l-4z 



2CIWa - 2 (2 C1C2 + N^Cl){Wx + W2/2) 



m 



-(1 - 2z)^{Ci + 2 C1C2 + iVeC2')g 



R 



(30) 



We again retain the dependence on quark masses in the above expression, including the terms 
proportional to m^. The operators in Eq. (30) are defined as 

QR=ipiSi)s+pipjSj)s+P , 

Wi^ — bi^"{l - 75)^aSi bj{l + 75)Sj , 

rub 



ml 



(31) 



where, as before, we have eliminated the color-rearranged operators Wj in favor of the operators 
Wi. The parametrization of the matrix elements of the above operators is given below, 

1 



{Bs\Qr\Bs)^-[2 



m., 1 



9 9 Ml 
fl Ml 



{Bs\Wi\Bs)=—il + 



nib V 2Nc 



mt 



{Bs\W2m = -i-l + — ] flM^ 



.9 9 I'Ml 

c2 T\/r2 B. 



I ^2 



1 "3, 



{BsmB,) = Ul + ^) flMl(^-l] a,, 



2Nr 



mt 
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(5,11^415.)=-- fsMl i^-^ - ij a,. (32) 

Note that in factorization approximation all the bag parameters ai should be set to 1. In addition 
to the set of kinetic corrections considered above, the effects of the interactions of the intermediate 
quarks with background gluon fields should also be included at this order. The contribution of 
those operators can be computed from the diagram of Fig. 2, resulting in 

+ 4 [P5 + Pe - ^7 - } ■ (33) 

The local four-quark operators in the above formulas are shown in Eq. (34): 

P,=kY{l - 75)5, 6,7^(1 - l^KG^^s, , 

Pg^i-^.^D'^^D V(l - 75)Sfe ttiG;, baail - 75)5, , 
P7=^&fcS'^^"7"(l - 75)«i tliGl, 6,7.(1 - 75)5^ , 

Ps-^bi'D''D''r{l - 75)Sfe tliGl, 6,7.(1 - 75)^,. (34) 

b 



m. 



Analogously to the previous section, and following Ref. [23], we parametrize the matrix elements 
in Eq. (34) as 



{B,\Pm = -flMl i^-^ - ij A. (35) 

We set j3i = 1 GeV^ to obtain a numerical estimate of this effect. It is clear that no precise 
prediction is possible with so many operators contributing to the lifetime difference. This, of 
course, is expected, as the number of contributing operators always increases significantly with each 
order in OPE. We can nonetheless evaluate the contribution of both l/m^ and 1/m^ by randomly 
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varying the parameters describing the matrix elements by ±30% around their "factorized" values. 
This way we obtain the interval of predictions of AF^^ and estimate the uncertainty of our result. 

C. Discussion 

Now we discuss the phenomenological implications of the results presented in the previous sec- 
tions. As usual in OPE-based calculations next-order corrections bring new unknown coefficients. 
In our numerical results we assume the value of the 6-quark pole mass to be nib — 4.8±0.2 GeV and 
/b^ = 230 ± 25 MeV. It might be advantageous to see what effects higher-order 1/ml corrections 
have on the value of AFb^. In order to see that we fix all perturbative parameters at the middle 
of their allowed ranges and show the dependence of AF^^ on non-perturbative parameters defined 
in Eqs. (28), (32), and (35): 

ATb, = [0.00055 + 0.17325, + 0.00245i - 0.0237^2 - 0.0024^3 - 0.0436^4 

+ 2 X lO-^ai + 4 X 10-^q;2 + 4 x lO'^a^ + 0.0009q;4 - 0.0007q;5 (36) 
+ 0.0002/3i - 0.0002/32 + 6 x lO'^ps - 6 x IQ-^P^ - 1 x lO^^/^g 
- 1 X 10-^/36 + 1 X 10-^/37 + 1 X 10-^ Ps] (ps~^). 

As one can see, 1/ml corrections provide rather minor overall impact on the calculation of AF^^. 
In particular, contributions of gluonic operators are essentially negligible. 

To obtain the complete Standard Model estimate of AF^^, we fix the perturbative scale in our 
calculations to fi = nih and vary the values of parameters of the matrix elements. Following [23] 
we adopt the statistical approach for presenting our results and generate 100000-point probability 
distributions of the lifetime, obtained by randomly varying our parameters within a ±30% interval 
around their "factorization" values. The decay constant /b^ and the b-quark pole mass nib are 
taken to vary within a la interval as indicated above. The results are presented in Fig. 3. This 
figure represents the main result of this paper [24] . 

There is no theoretically-consistent way to translate the histogram of Figure 3 into numerical 
predictions for AF^^ • As a useful estimate we give a numerical prediction by estimating the width 
of the distribution Fig. 3 at the middle of its height and position of the maximum of the curve as 
the most probable value. We caution that predictions obtained this way should be treated with 
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FIG. 3: Histogram showing the random distribution around 
the central values of various parameters of Eqs. (24, 26, 30, 
33) contributing to Bg-lifetime difference AFb^. 

care, as it is not expected that the theoretical predictions are distributed according to the Gaussian 
distribution. Nevertheless, following the procedure described above one obtains 



AFs, = 0.072iHl^ ps-\ 



0.104 ±0.049, 



(37) 



where we added the experimental error from the determination of Tg and theoretical error from 
our calculation of AF^^ in quadrature. 



IV. NEW PHYSICS CONTRIBUTIONS TO LIFETIME DIFFERENCE 

In the previous section we have shown that 1/m^-corrections to the lifetime difference of the 
light and heavy cigenstates in the Bg system are quite small, which makes the prediction of AF^^ 
quite reliable ^. Additionally improving the accuracy of the lattice or QCD sum rule determinations 
of non-perturbative "bag parameters" in Eq. (36) would make this prediction even more solid. 

^ As was argued in Ref. [16], perturbative scale dependence can be further reduced by switching to a different basis 
of leading-order operators. 
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In this respect, it might be interesting to consider the effects of New Physics on the hfetime 
difference in Bg system. Why would it be worthwhile to perform this exercise, especially since it is 
known that AFs, is dominated by the on-shell, real intermediate states? Wouldn't AB — 1 New 
Physics amplitudes that can potentially affect AP^^ already show up in the experimental studies 
of exclusive Bs decays? This is indeed so. However, it might be difficult to separate New Physics 
effects from the dominant (but somewhat uncertain) Standard Model contributions, as theoretical 
control over soft QCD effects is harder to achieve in the calculations of exclusive decays despite 
recent significant advances in this area [25]. 

It was recently pointed out that NP contributions can dominate lifetime difference in — 
system in the flavor SU{3) limit [26]. In that system this effect can be traced to the fact that 
the SM contribution vanishes in that limit. While similar effect does not occur in Bg mixing, 
good theoretical control over non-perturbative uncertainties in the calculation of AP^^ makes 
calculations of NP contributions worthwhile. In S^-system one can show that 

AMb^ = 2|Mi2|, 

^^^^ - amb. ■ ^^^^ 

In the Standard Model the phase difference between the mixing amplitude and the dominant decay 
amplitudes is arg {—V*Jycs/V^^Vts), i.e. essentially zero. If NP contribution has a CP- violating phase 
that exceeds that of the Standard Model, one can write, denoting 2^ = arg(Mi2Pi2), 

APb^ = 2|Pi2|cos2^. (39) 

Since in the Standard Model ri2 is dominated by the h — > ccs transition, its phase is negligible. 
Then, as was pointed out long time ago [19, 20], CP- violating contributions to M12 must reduce 
the lifetime difference in S^-system, 

APs^^APlf cos2e, (40) 

where 2^ is a CP- violating phase of M12, which is assumed to be dominated by some AB — 2 New 
Physics. 

Contrary to CP- violating AB = 2 NP contributions to M12, any AB = 1 NP amplitudes 
can interfere with the Standard Model ones both constructively and destructively, depending on 
the model. Since no spectacular NP phases have been observed in Bg mixing, it appears that 
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Mi2 is dominated by the Standard Model CP-conserving contribution. In that case, the phase 
arg (Mi2r^2) = arg (r^2) = 2^' is dominated by the phase of New Physics contribution to r^2. In 
that case 

Ar5^ = Ar|f + Ar^fcos2e', (41) 

where AF^^^ is a contribution resulting form the interference of the SM and NP AS = 1 operators, 
which can either enhance or suppress AF^^ compared to the Standard Model contribution. We 
shall compute AF^^^ by first employing the generic set of effective operators, and then specifying 

to particular extensions of the SM. We shall concentrate on CP-conserving contributions. 

Using the completeness relation the NP contribution to the -8°-^ lifetime difference becomes 



^ {Bs\lmT\Bs) , where (42) 



where we represent the generic NP AB — 1 Hamiltonian Ti^f as 

n^p=' = E IWQi + C2{fj^)Q2\ , (43) 



where the spin matrices Fi 2 can have an arbitrary Dirac structure, Dqqi are some New Physics- 
generated coefficient functions [26], and Ci 2(Ai) are Wilson coefficients evaluated at the energy 
scale ^l. This gives us the following contribution to the lifetime difference: 
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^^b! = ^Y.D,,'V:,V,,, {K,5,M + K^5,,5u) E h{x,x'){B,\Of'\Bs), (44) 



where k, I are the color indices, {Ka} are combinations of Wilson coefficients, 

Ki = {C2C2N, + (C2C1 + C2C1) ) , K2 = CiCi (45) 

with the number of colors Nc — 3, and the operators O^'^' are the following: 

Of'^(hr^Yr2Si) (6,Fi7,F.s,) 
of =(6,F- ^T2Si) (6fcFi i>F.6,) 

Of =(6,F-F2S,) (6,Fi ^r^sj) , (46) 
14 



where ^ is the 6-quark momentum operator. Defining Zg = mi /ml and Zqi = mi, /ml the coeffi- 



cients Ij{zq,Zq') can be written as follows: 



k*mb 



h{Zq,Zg>)-- 



487r 
k* 



24mb7r 



1 - 2{Zq + Zq') + {Zq - Zq'f 
l + {Zq + Zq') - 2{Zq - Zq'f 



k* 



h{Zq, Zq')^ — ^/z^ [l + Zq' - Zq] , 

k* 



h{Zq, Zq')^- — y/z^ [1 - Zq' + 



(47) 



h{Zq, Zq'] 



Stt 

k*mi, 



\J ZqZqi , 



where k* — (mb/2) [1 — 2{zq + Zqi) + {zq — Zq')'^]^^'^ . This is the most general formula for the New 
Physics contribution to the lifetime difference in Bg mesons. We now look into two particular 
examples extensions of the Standard Model, multi-doublet Higgs models and Left-Right Symmetric 
Models, that can contribute to AF^^. 



A. Multi-Higgs model 

One of possible realizations of New Physics is a multi- Higgs doublet model [27]. Many of SM 
extensions, particularly the supersymmetric ones, require extended Higgs sector in order to break 
additional symmetries of NP down to SU{2)l x U{1) of the Standard Model. These constructions 
contain charged Higgs bosons as parts of the extended Higgs sector. These models provide new 
flavor-changing interactions mediated by charged Higgs bosons, which lead to rich low-energy 
phenomenology [28, 29]. In the low-energy limit, charged Higgs exchange leads to the following 
four-fermion interaction [30], 

-H^hh' = b^T^ql q^T.s, , (48) 

where Fj, i = 1, 2, are 

Ti^mbV^cotpPL - mcV^tenipPji, 
15 



r2^msVcs cot PPr - rricVcs tan ^Pl, (49) 

where Pl,r = (1 T 72)/2- Inserting Eq.(48) into Eq. (44) leads to a contribution to the hfetime 
difference {ArBj^s)chH from three operators with various coefficients, 



F 6 I c6 cs) j^^^^ (4X2^/1 cot^ P + 2{coe (5ml^sh - m,^,h){K2 - K^)) 
+ {Q2) (-2/^1^/1 cot' /? + (cot' Pml^j2 - rrib^ch) {K2 - 



ChH MbTb. Ml 



+ (Qs) (^1 + K2) (zc tan' (5h - . (50) 

7j = Ii{zc.i Zc), Ki are defined above, and (Qi) are 

For values of Mh = 85 Gel/ and cot (3 = 0.05 [11] we obtain (Ar^^r,)^^^ « 0.006. This is about 
6% of the Standard Model value, too small to constrain the model from this observable. The 
dependence of {ArBj^s)chH mass of the Higgs boson is given in Fig. 4. 



B. Left-Right Symmetric Models 

One of the puzzling features of the Standard Model is the left-handed structure of the elec- 
troweak interactions. A possible extension of the SM, a Left-Right Symmetric Model (LRSM) 
assumes the extended SU{2)l x SU{2)ji symmetry of the theory, which restores parity at high 
energies [31]. While in the simplest reahzations of LRSM the right-handed symmetry is broken at 
a very high scale, models can be consistently modified to yield Wfl-bosons whose masses are not 
far above 1 TeV range [32]. In this case fiavor-changing interaction from l^i^-bosons can affect 
ATb^ (for a similar effect in D-mixing, see [26]). 

In principle, manifest left-right symmetry requires that couplings to left-handed particles to 
be the same as the once to the right-handed ones, e.g. ql — Qr. This also assumes that the 
right-handed CKM matrix V^^^ should be the same as the left-handed CKM matrix Vi^. In this 
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FIG. 4: Dependence of ychH on the mass of the Higgs boson. 
SoUd Une: tan/3 = 20; dashed Une: tan/3 = 10; dotted hne: 
tan/3 = 5; dash-dotted hne: tan/3 = 3. 



case, kaon mixing constraints exclude My/n < 1-6 TeV [33] (direct constraints are weaker by 
approximately factor if two). However, v}j^^ could also be quite different from the Vik, as long 
as it is still unitary. In this case of non-manifest left-right symmetry the bounds on are 
significantly weaker, Mw^ > 0.3 TeV from kaon mixing [34]. To assess the contribution from Wr 
to Ar^^, we equate 



(R) 



^qq' — ^cb ^cs ^2 



cb 



ri,2 = i^Pi 



R 



(52) 



in Eq. (44) and evaluate the respective operators. Here / = gj^/SM^P'^ , and we assume 
Qn = kql. In the studies of non-manifest LRSM, we shall also assume k = 1, 1.5,2 [35]. At the 
end, LRSM gives the following contribution to the value of ATbJTb/- 

2 



LR 



w 



(R) 
W ■ 



\Cx{Q2)-'^C^m\- (53) 



The dependence of (Ar^^/r^Jj^/j on the mass of the Wn boson is given in Fig. 5. We see that 
contrary to the D-meson case [26, 36], i3s-mixing could provide decent constraints on the values of 



Mp^ . For instance, in a non-manifest LRSM (with relevant V^j 
one obtains i^IS^ bJ'^ b^lr — —0.04 This is a rather large contribution to AF^^, more than a third 
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(a)Dependence of ATbJ^s on Myy ■ Solid line: 



non- manifest LRSM (k = 1), dashed line: 
manifest LRSM. 



(b)Dependence of AFs, /Fs on M^^' in 



non-manifest LRSM. Solid line: k — 1, dashed 
line: n = 1.5, dotted line: k — 2. 

FIG. 5: Contributions to ATb^/^s in the Left-Right Symmetric Models. 



of the absolute value of the Standard Model contribution and of the opposite sign. The LRSM 



contributions for k > 1 are even larger. As expected, in the case of manifest LRSM = Vi 

the contribution from this model is less marked, {ATbJTbJlr. < 0.002 for M^y > 800 GeV 



V. CONCLUSIONS 

We computed the subleading 1/ml corrections to the difference in the lifetimes of Bg mesons. 
We showed that they can be parameterized by 13 nonperturbative parameters, which we denote a, 
and Pi. We adopted the statistical approach for presenting our results and generate 100000-point 
probability distributions of the lifetime difference, obtained by randomly varying our parameters 
within a ±30% interval around their "factorization" values, except for the case when the parameters 
are known from lattice QCD. In this case they are taken to vary within a la interval as indicated 
above. 

The results are presented in Fig. (3). While there is no theoretically-consistent way to translate 
the histogram of Fig. 3 into numerical predictions for ATbJTs, we provide an estimate by taking 
the width of the distribution Fig. 3 at the middle of its height as 1-cr variance and position of the 
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maximum of the curve as the most probable value, 

Ar 

ATb, = 0.0721^30 Ps"\ -jr^ = 0-104 ± 0.049, (54) 

The effects of 1/m^ corrections to calculations of AF^^ are shown to be small. 

We also looked into AB = 1 New Physics contribution to the lifetime difference in the Bg 
system. We have shown that these contributions can both enhance or reduce the Standard Model 
contribution. We considered the most general four-fermion effective Hamiltonian, which can be 
generated by any reasonable extension of the Standard Model and derived its contribution to AF^^. 
We then evaluated effects of charged Higgses and right-handed Ws on the lifetime difference. While 
the contribution of charged Higgs was shown to be negligible in AF^^, LRSM can be constrained 
with measurement of AF , provided lattice or QCD sum rule community provide better estimates 
of non-perturbative parameters entering the SM calculation of the lifetime difference in Bg mesons. 
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